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C*") ' We study a variational model for a diblock copolymer- homopolymer blend. The energy 

functional is a sharp-interface limit of a generalisation of the Ohta-Kawasaki energy. In one 
dimension, on the real line and on the torus, we prove existence of minimisers of this functional 
and we describe in complete detail the structure and energy of stationary points. Furthermore 
we characterise the conditions under which the minimisers may be non-unique. 
*I3 ' In higher dimensions we construct lower and upper bounds on the energy of minimisers, 

and explicitly compute the energy of spherically symmetric configurations. 
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■ 1 Introduction 

■ 1.1 Micro-phase separation 

' In this paper we study the functional 

o 



co/ \V(u + v)\+c u \\7u\+c v \Vv\ + \\u-v\\jj- imN) if (u,v) € K x , 
Fi(u,v) = { Jm n Jm n Jwl n 

oo otherwise, 



h ■ (i) 



where the coefficients Ci are nonnegative (not all equal to zero) and0 

K x := \ (u,v) € (BVQR^)) 2 : u(x),v(x) £ {0, 1} a.c., and uv = a.e., and 



v 



Under the additional constraint u + v = 1, this functional is the sharp- interface limit of a well- 
studied variat i onal model for melts of diblock copolym ers [ChoOll ICR031 [CS061 IFH011 IMur021 
IRW00I IRW02I IRWOSal IRWOSbl IRW051 IRWOBal IRWOBbj. T his underlying diffuse interface model 
is also closely related to the functional studied in |Miil93j . Such polymers consist of two parts, 
labelled the U and V parts, whose volume fractions are represented by the variables u and v. 
The U and V parts of the polymers repel each other, and this repulsion leads to micro-phase 
separation: phase separation at a length scale comparable to the length of a single molecule. The 
case studied here is known as the strong segregation limit, |BF99j , in which strong repulsion causes 
strong demixing of the constituents — hence the restriction of K x to characteristic functions. The 
modeling assumption here is that stationary points of F x under constrained (i.e. fixed) mass J RN u, 
in particular minimisers, represent the structures formed by the polymers. 



^^Where we do not explicitly specify the integration measure, we use the Lebesgue measure. 
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Although the various simplifications leading to F\ have obscured the connection between this 
functional and single molecules, the character of the various terms is still recognisable. The 
intcrfacial penalisation terms, i.e. the first three terms, are what remains of the repulsion in the 
strong segregation limit, and these terms favour large-scale demixing. The last term \\u— v\\ H -i. 
on the other hand, penalises such large-scale separation and arises from the chemical bond between 
the U and V parts of the polymer molecules. 

These competing tendencies cause the functional F\ to prefer structures with a specific length 
scale, as we now illustrate with a simple example in one space dimension. For simplicity we take 
as spatial domain the unit torus Ti, i.e. the set [0, 1] with periodic boundary conditions; all global 
minimisers under the condition u + v = 1 then are of the form shown in Figure Hal 
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(a) n-block structure on the unit torus Ti , with u+v = 1 



(b) n-block structure on K or Ti 



Figure 1: 



For such structures the value of the functional is 

Ft = 2n{c u + c v ) + — ^ 

as can be seen from the results in Section [3J If we consider c u and c v to be fixed, the energy F\ 
is clearly minimised at a finite value of n. When we study the one-dimensional case on R without 
the restriction u + v = 1 in more detail, in Section [3l we shall see that the energy actually favours 
a specific block width rather than a specific number of blocks. 



1.2 Blends of co- and homopolymers 

For u + v ^ 1, F\ is a model for blends, mixtures of diblock copolymers and homopolymers; the 
homopolymer is considered to fill the space not occupied by the diblock copolymer and has local 
volume fraction 1 — u — v. 

The inclusion of homopolymers into a block copolymer melt opens the possibility of structures 
with two distinct length scales. The repulsion between the two blocks creates micro-phase sepa- 
ration at the length scale of the polymer, as described above. At a larger length scale structures 
are observed in which regions of pure homopolymer and pure copolymer alternate. 

Blend systems show a tremendous wealth of behaviour. For instance, many different types 
of macrodomain geometry have been obse rved: sphe res [KHH94| ION97[ IUD05| IZJM05j . cylin- 
ders [KWT88], dum bbells [01 95], helices |HMY+01] . labyrinths and sponges [LSHG941 IPcoM 
OI95J, ball-of-thrcad TSHG94 , and many more. In addition, the microdomains have varying ori- 
entation with respect to this macrodomain geometry. In many cases the micro- and macrodomain 
geometry appear to be coupled in ways that are not yet understood. 

There is extensive literature on such blend systems, which is mostly experimental or numerical. 
For the numerical experiments it is de rigeur to apply a self-consistent mean field theory and 
obtain a generalisation of the Ohta-Kawasaki |OK86j mode l (see e .g. |NH831 ION971 ICR05j ). Of 
the resulting model the energy F\ is a sharp- interface limit |Bal901 ICR05] . 

At the level of mathematical analysis, however, little is known. What form do global and local 
minimisers of F\ take? (Do they even exist? The issue of existence of global minimisers of F\ 
on K is first addressed in this paper.) Docs the functional indeed have a preference for layered 
structures, as the numerical experiments suggest? What structure and form can macrodomains 
have? Can we observe in this simplified functional F\ the breadth of behaviour that is observed 
in experiments? All these questions are open, and in this paper we provide some first answers. 
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1.3 Results: global minimisers in one dimension under constrained mass 

The first part of the paper focuses on the one-dimensional situation. 

1.3.1 Existence 

The existence of global minimisers under the constraint of fixed mass follows mostly from classical 
arguments (proof of Theorem 13. 7[) . The non-compactness of the set R can be remedied with the 
cut-and-paste techniques that we introduce to study non-uniqueness (sec below). 

One non-trivial issue arises when e.g. cq — c u = 0, in which case the functional F\ provides no 
control on the regularity of u. We obtain weak convergence in L 2 for a minimising sequence, and 
therefore a priori we can only conclude that the value set of the limit functions is [0, 1], the convex 
hull of {0, 1}; as a result the limit (u,v) need not be an element of K\. With a detailed study of 
the stationarity conditions on u we show that stationary points of F\ only assume the extremal 
values and 1. The existence of a minimiser then follows from standard lower semi-continuity 
arguments. 

1.3.2 Characterisation of macrodomains 

In the one-dimensional situation a macrodomain is a finite sequence of alternating U- and V- 
'blocks' or 'layers' as in Figure [TBI Choksi and Ren [CR05] studied such macrodomains defined on 
the torus of length L, but their techniques apply unchanged to the real line also. They showed 
that if such a macrodomain is stationary, then all interior blocks have equal width, while the end 
blocks are thinner. 

The exact dimensions of the blocks arc fully determined by the number of blocks, the total 
mass, and, in the case of the torus, the size of the domain (see Theorems 13.21 and 13. 3p . It is 
instructive to minimise F\ within classes defined by a specific choice of the sequence of U- and 
V-blocks; Figure [2] shows this minimal energy for different classes and different values of the mass. 

El 
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2 4 6 8 

Figure 2: Energy per unit mass for the one-dimensional case R, according to the calculations in 
Section l3~5l M is the total U-mass; for the surface tension parameters (see Lemma |2~3| the values 
d u o = ljrf«u = 0.7 and d V Q = 0.3 are chosen. The graphs belong to the following structures, as 
indicated in the figure as well (the lighter coloured blocks are V-blocks, the darker ones U-blocks): 
(a) UV and VU, (b) UVU, (c) VUV, (d) UVUV and VUVU, (e) UVUVU, (f) VUVUV. The circle 
indicates where the optimal structure changes. 



3 



1.3.3 Characterisation of constrained minimisers 

We extend the results of Choksi & Ren into a full characterisation of global minimisers, by showing 
that there exists a global minimiser with only one macrodomain, and by fully characterising all 
other global minimisers in terms of the parameters and the morphology (Theorem 13. 7)) . 

This characterisation shows that non-uniqueness of minimisers can take two different forms. 
The first is the possibility that two different UV-sequences with the same mass have the same 
energy, as is illustrated by the encircled intersection in Figure O This is a common occurrence in 
variational problems, where a parameter change causes the global minimum to switch from one 
local minimiser to another. 

The second form of non- uniqueness is related to the fact, which we prove in Section [3.3[ that 
two separate macrodomains can be translated towards each other and joined together without 
increasing the energy. In fact, in many cases the energy strictly decreases, and it is this possibility 
of strict decrease that allows us to rule out many cases. This leaves us with a set of conditions for 
the case of unchanged energy that must be fulfilled when a non-unique global minimiser contains 
more than one macrodomain (see Theorem 13 - 7[) . This type of non-uniqueness is specific for the 
problem at hand, and produces not a discrete set of minimisers but a continuum, parametrised by 
the spacing between the macrodomains. 

Although the focus of this paper lies on the unbounded domains K and R JV , we make a brief 
excursion to extend the characterisation of global minimisers to the case of the torus Tl with 
length L (Theorem EH)]). 



1.3.4 A lower bound 

Figure [2] and more clearly Figure [3] illustrate that as the imposed mass increases the number of 
blocks of the global minimiser(s) also increases. In Section [331 wc calculate values of the energy 
for various global minimisers, and show that the thickness of the internal layers approaches the 
optimal spacing of 

2m -.^^{cu + c^l^ 

for M — > oo while the width of the end layers converges to half this value (Remark 13. 12)) . 

As a corollary we obtain an explicit and sharp lower bound for the energy on R (Theorem l3.11[) : 

1/3 

(2) 



Fi(u,v) > 2(c + mm(c u ,c v )) + ( - j (c u + c v fi 3 



The fact that the lower bound is sharp is significant. For instance, the affine dependence of the 
right-hand side on the mass J u implies that the minimal energy per unit of mass, F\ (u, v) (J u) , 
is generically not attained at any finite mass, but only in the limit J u — > oo. 

The word 'generic' refers here to the assumption that cq + min(c u , c v ) > 0, and the alternative 
case Co = c u — (or cq = c v = 0) is fundamentally different. In this latter case macrodomains can 
be split and joined without changing the energy. 

The characterisation of global minimisers also allows us to establish an asymptotically sharp 
upper bound (Theorem 13. 1 1[) : 

Jim^inf ^M^F^v) : (u,v) e K u J u = m| = (^j ' {c u + c v ) 2 '\ (3) 
In the limit M — > oo, the bound @ coincides with ([5]). 



1.4 Results: higher dimensions 
1.4.1 Energy bounds 

A common strategy in the study of pattern-forming systems is not to make any Ansatze about 
the morphology but to search for weaker characterisations of behaviour. As an example of this 
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in the field of block copolymers, Choksi proves that for pure diblock melts the energy is bounded 
from below by a lower bound with a certain scaling in terms of the physical parameters — without 
making any a priori assumptions on the morphology |Cho01| . This scaling is shared by periodic 
lamellar structures with a specific lamellar separation. 

For the case at hand, the one-dimensional analysis provides both a lower and an upper bound 
on the energy in one dimension. Weakening the lower bound (J2J to 

Fi{u t v)> ^ (c u + c v ) 2 / 3 J u, (4) 

one might conjecture that the lower bound ^ holds in R , again without making any a priori 
assumption on the morphology. However, we have no proof of this conjecture, and in fact, results 
on mono- and bilaycr stability (see Section 11.4.21 below) suggest that such a conjecture may only 
hold for certain choices of the parameters. In Section 0] we instead prove a lower bound which is 
also linear in mass, but has a smaller constant (Thcorcm l4.ip . 

The explicit construction used to prove the upper bound |3} suggests a natural strategy for 
proving a similar upper bound in M. N . In Sections 14.21 and 14.31 we extend one-dimensional min- 
imiscrs as lamellar structures in M. , and prove the same upper bound ((3|) in M. N . Here the main 
step in the proof is the estimation that 'boundary effects' as a result of the cutoff to finite mass 
are of lower order. 

In Section 14.31 we use the same idea to calculate the energy values of some structures with 
spherical geometry: either solid spheres of one phase (U or V) surrounded by a spherical layer of 
the other phase (micelles), or ring-shaped layered structures. In both cases the asymptotic energy 
exceeds that given by the upper bound (J3j) , indicating that they can not be global minimisers 
in R N . 



1.4.2 Monolayer and bilayer stability in periodic strips 

In a companion paper |GP07j we study the stability with respect to a certain class of perturbations 
of monolayers and bilayers, i.e. straight layered structures with one respectively two lines of U-V 
interface, in a periodic strip Tl x K. There we show that for sufficiently large L a monolayer (the 
simplest lamellar structure, of the form UV) is always unstable, while the stability of a bilayer 
(UVU or VUV) depends on the parameters. For the case of a UVU bilayer with optimal thickness, 
for instance, we prove a stability criterion of the form 

, c u + c v f L 

stability > g 



c + 2c u + c v \ (c + 2c„ + Cu) 1 / 3 1 

where g is a continuous function with values in (0, 1). Therefore, the bilayer can be stable or 
unstable, depending on the relative values in the interface penalisation parameters. Note that the 
relative value of c u + c v should not be too small in order to have stability. More about the special 
role of c u + c v follows in Section 12.21 



1.5 Related work: partial localisation 

In previous work, one of the authors (Pcletier) and Roger studied a related functional whose deriva- 
tion was inspired by lipid bilayers [PR06] . Lipid bilayers might be considered block copolymers, 
and therefore it is not surprising that the functional considered in PR06J is similar to F\ : 

!e I |Vm| + -di(u, v) if (u,v) e IC £ , 
Jr 2 £ (5) 

oo otherwise. 

Here u is the volume fraction or density of lipid heads, v is the volume fraction of lipid tails, g?i(-, •) 
is the Monge-Kantorovich distance and 

K e := \ (u,v) G BV(M 2 ;{0, 1/e}) 2 : uv = a.c., and u= v = M \ . 
I Jr 2 Jr 2 J 
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Apart from the choices Co = c v = and c u = 1, the main difference between ([T]) and ([S]) is the 
different non-local term. 

Note that the scaling (constant mass but increasing amplitude l/e) implies that the supports 
of u and v shrink to zero measure. The main goal in |PR06j was to investigate the limit e — * 0, 
and connect the limit behaviour to macroscopic mechanical properties of the lipid bilaycrs such 
as stretching, bending, and fracture. 

The authors studied sequences (u £ ,v £ ) for which the rescaled energy Q £ := e~ 2 (T £ — 2M) 
remains finite. They revealed a remarkable property of the functional Q £ (or T £ ): boundedness of 
Q £ {u £ , v £ ) implies that the support of u £ and v £ becomes close, in the sense of Hausdorff distance 
between sets, to a collection of closed curves of total length M . The curve-like behaviour indicates 
partial localisation: localisation in one direction (normal to the limit curve) and non-localisation 
in the direction of the tangent. In addition one can recognise resistance to stretching (because of 
the fixed length) and resistance to fracture (because the curves are closed). Moreover, the curves' 
support is approximately of 'thickness' 2s, indicating an underlying bilayer structure. The authors 
also showed that Q £ Gamma-converges to the Elastica functional, which penalises the curvature 
of curves, showing a tendency of the limit curves to resist bending. 

These results suggest considering similar limits for the functional F\. In fact the subscript 1 
in F\ and K\ already refers to the appropriate rescaling: 



F £ {u,v) 



el cq I \V(u + v)\ + c u I \Vu\+c v I |V«| ) + -||u - v\\%^ if (u, v) G K e , 



oo otherwise, 

(6) 

where 

K £ := < (u, v) G (BV(M Ar )) 2 : u(x), v(x) G {0, l/e} a.e., and uv — a.e., and / u = I v 
I Js. N Jr n 

As mentioned above, in the companion paper |GP07j we investigate the stability of bilayers, and 
show that parameter choices exist for which they are stable: this provides another suggestion that 
the functional F £ may display similar behaviour in the limit e J, 0. This is work for future research. 



2 Preliminary definitions 

2.1 Problem setting 

In this paper we mostly consider as domain the whole space R w ; however, sometimes we will make 
an excursion to the torus , i.e. a periodic cell rifc=i[0i^i] with the endpoints of each interval 
identified. 

Definition 2.1. For f 6 L 1 (R Ar ) (or L 1 (T^ r ) y ) with J f = and compact support, 

11/11*-! := JfG*f, (7) 

where G is a Green's function of the operator —A on W N (or on T^). We define the space 
H^ 1 (M. N ) as the completion of {/ G L^(M. N ) : supp / compact, J RN f = 0} with respect to the norm 

in Similarly i? _1 (T^ r ) is defined as the completion of |/ G L 1 (T^ r ) : J jN f = o| with respect 
to this norm. 

On the zero average condition of / is necessary in order for G * f to respect the topology 
of the torus: 

/"/ = -/ A(G * /) = 0. 
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This condition also allows for a convenient reformulation of the norm ([7]) in terms of the Poisson 
potential (f>t of /, given by 

cj> f = G*f, 

such that 

\\f\\ 2 H-i= ff<t>f = I W/1 2 . (8) 



Fi(u,v) 



In some cases it will be useful to add a constant to <f>f\ note that this can be done without changing 
the value in ([5]). 

If the set Hq is defined as the completion of C^.(R N ) (or C 1 (T^ r ) with zero mean) with respect 
to the norm ||g||^i = J |Vg| 2 , then (J7J) is the dual norm of Hq with respect to the i 2 -inner product 
and satisfies 

J fg<\\f\\H-4g\\Hi, 

for all / e H~ x and g e Hq. 

We repeat the definition of F\ and K\ for convenience. 

Definition 2.2. Let cq, c U) and c v be real numbers. Define 

c \V{u + v)\+c u \Vu\ + c v |Vt>| +||«-t;||| r _, if {u,v) £ K x , 

oo otherwise, 
where the admissible set is given by 

) e (BV{R N )) : u{x), v(x) £ {0, 1} a.e., and uv = a.e., ana 

We will require Co, c u and c v to be non-negative and assume that at least one of these coefficients 
is positive. See also Remark 12.51 

Sometimes we consider the case of the torus instead of E . It is understood that in the above 
definition the instances of arc then replaced by . 

Another, equivalent, form of the functional will be useful, in which the penalisation of the three 
types of interface U-0, V-0, and U-V, is given explicitly by surface tension coefficients dki- 

Lemma 2.3. Let the surface tension coefficients be given by 



d u o 
dim 



— Cu + CO, 
= Cv + C , 
Cu l C v , 



N on- negativity of the Ci is equivalent to the conditions^ 

< dki < dkj + dji for each k =/= I. (9) 

Then 

F (uv) = { d ^ nN ^( S ^) + dvoH"- 1 ^) + dvvH"- 1 ^) + \\u- v\\ 2 H ^ if (u,v) g Ki, 
1 1 oo otherwise. 

where Ski is the interface between the phases k and I: 

S u o = d* supp u\d* supp v, 
Svo = d* supp v\d* supp u, 
S U v = d* supp u n d* supp v, 

and d* is the essential boundary of a set. 



2 The indices j, k,l take values in {u,v,0} and the d^i are taken symmetric in their indices, i.e. d vu := d uv etc. 
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Remark 2.4. The essential boundary of a set consists of all points in the set that have a density 
other than or 1 in the set. Details can be found in [AFP001 Chapter 3.5]. 



Proof of Lemma \2.3l The main step in recognising the equivalence of both forms of F\ is noticing 
that, if u is a characteristic function, then 



Note the different interpretations of the coefficients Cj and the surface tension coefficients dki- 
The latter have a direct physical interpretation: they determine the mutual repulsion between the 
different constituents of the diblock copolymer-homopolymer blend. For example, the value of d uv 
(as compared to the values of d U Q, d v o and 1, the coefficient in front of the i7 _1 -norm) determines 
the energy penalty associated with close proximity of U- and V-polymcrs. In particular, if one 
of these surface tension coefficients is zero, the corresponding polymers do not repel each other 
and many interfaces between their respective phases in the model can be expected. On the other 
hand the coefficients Ci, when taken separately, do not convey complete information about the 
penalisation of the boundary of a phase. If for instance c u = 0, but c v ^ 0, the part of the 
U-phase interface that borders on the V-phase still receives a penalty, because d uv = c v . For this 
reason the use of surface tension coefficients makes more sense from a physical point of view. For 
the mathematics it is often easier to use the formulation in terms of Ci . 

Remark 2.5. The condition (HJ) can be understood in several ways. If, for instance, d uv > 
d u o + d v o, then the U-V type interface, which is penalised with a weight of d uv , is unstable, for 
the energy can be reduced by slightly separating the U and V regions and creating a thin zone 
of inbetween. A different way of seeing the necessity of ([9]) is by remarking that the equivalent 
requirement of non-negativity of the Ci is necessary for F\ to be lower scmicontinuous in e.g. the 
L 1 topology. Our assumption that at least one Ci is positive is equivalent to assuming that at least 
two dki are positive. 



The behaviour of the model described by F% is crucially different in the two cases d uv > (c u +c v > 
0) and d uv — (c u = c v = 0). The statements made in the introduction such as 'the functional 
Fi prefers structures with a definite length scale' actually only hold in the case d uv > 0. For most 
results in this work we will assume this condition to hold, and to justify this we now show with 
an example how the case d uv = is different. 

Consider the one-dimensional case, take f2 to be the torus Ti, and fix Co = 1 and c u — c v = 0, 
or equivalcntly d uv = and d u o = d v o — 1. Restricting ourselves to functions (u, v) £ K\ with 
Jq 1 u = Jq v = M, for some fixed mass < M < 1/2, we find that for any (u, v) there are at least 
two U-0 or V-0 type transitions, and therefore 



On the other hand, equality is only reached if u — v = 0, which is not possible for positive mass M . 
But the value 2 can be reached by a sequence of approximating pairs (u n ,v n ), 




□ 



2.2 The role of d, 



'UV 





\x\ < n and 
otherwise 



\x\ < n and — < x < 2khi f or some k <E Z 
otherwise 



< x < 2£ for some k E Z 



Then (u n ,v n ) € K\ and 
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• fo\(un + v n y\ = f 1 \ x \_ n>n] \ = 2; 



• In Section [3] it is calculated that a single one-dimensional monolayer of width 2m and height 
1 satisfies | j it — = 2m 3 /3; extending this result to the functions (u n ,v n ), which are 

concatenations of n 2 such monolayers, each of width 2/n, we find ||tt — = n 2 -2n~ 3 /3 = 



Consequently, Fi(u n ,v n ) converges to 2 for n — > oo. 

This sequence illustrates the preferred behaviour when d uv = 0: since the interfaces between 
the U- and V-phases are not penalised, rapid alternation of U- and V-phasc effectively eliminates 
the H -norm, reducing the energy to the interfacial energy associated with a single field u + v. 

3 Global minimisers in one dimension 

In this section we fully characterise the set of global minimisers of Fi in one space dimension, i.e. 
N = 1. Our main discussion concerns the case of R, but in Section f3T4l we will briefly mention 
results on the torus Tl- 

In one space dimension it is useful to regard admissible functions (u, v) as a sequence of 
blocks. A U-block, a V -block, and a 0-block are connected components of suppu, suppu, and 
R \ supp(u + v), respectively. Adjacent blocks are separated by transitions or interfaces. We will 
sec below (Corollary 13. 6| that any stationary point has a finite number of interfaces, even if either 
d U Q or d V Q vanishes. 

If (u, v) is an admissible pair, each of the connected components of its support supp(u + v) is 
in fact a macrodomain in the sense of the introduction. If there is only one such macrodomain, we 
call the configuration connected. Thinking about the structures in terms of sequences of blocks, 
we can specify connected configurations up to block width and translation by a sequence of U's 
and V's, e.g. UVUVU. 

Characterising the set of global minimisers falls apart into two steps: 

A For a given macrodomain we describe the optimal spacing between the transitions; 

B We derive necessary conditions for the occurrence of a disconnected global minimiser, i.e. a 
global minimiser with more than one macrodomain. 

In addition we use the techniques of part B above to prove the existence of a global minimiser. 

In Section [3~2l we first describe the characterisation given by Choksi and Ren |CR05j of the in- 
ternal structure of macrodomains, which essentially coincides with part A above. We then continue 
in Section [3.31 by showing that the support can be reduced to a single connected component; this 
also provides necessary and sufficient conditions for non-uniqueness (Theorem 13. 7[) . The reduction 
to a single macrodomain also allows us to prove an existence result (Theorem 13. 7[) . Finally, in 
Section 13.51 we calculate the values of these minimisers and derive a lower bound for the energy 
per unit of mass. 

3.1 Stationarity 

Because the set of admissible functions K\ is not locally convex we need to carefully formulate 
the notion of stationary point. 

Definition 3.1. We call (u,v) € K\ a stationary point of F\ if for any sequence (u n ,v n ) C Ki 
such that u n —> u in L 1 and !)„-m in L 1 , 



27J" 1 
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= 0. 

t=o 



As a consequence of this definition, if t i— > (it(i), i>(i)) is a curve in Xi, with (u(0),v(Q)) a 
stationary point of F\ , then 

In the proofs of the results in Section [3~2l a special case of this is used: for a connected configuration 
in one dimension that is stationary under constrained mass the derivative of F\ with respect to 
mass-preserving changes in the position of the interfaces is zero. 

3.2 Characterisation of macrodomains 

For periodic domains, Choksi and Ren [CR05] have given a characterisation of the structure of 
macrodomains. For its formulation it is useful to define three types of interface. Interfaces 0-U 
and U-0 interfaces are considered to be of the same type, as are 0-V and V-0 interfaces and U-V 
and V-U interfaces. Choksi and Ren's conclusions are 

Theorem 3.2 ( [CR0 5]). Let (u,v) be a stationary point of Fi on the torus Tl under constrained 
mass, with supp(u + v) connected and with a finite number of interfaces. Then 

1. Each pair of adjacent U-V type transitions is separated by the same amount; i.e. each U- or 
V -block is of the same width, with the exception of the two end blocks. 

2. In the cases UVUV. . . U and VUVU. . . V the end blocks are half as wide as the internal 
blocks. 

3. In the case UVUV. . . V (or the mirrored configuration VUVU. . . U) there is an additional 
relation that determines the width of the end blocks. 

The case of R was not explicitly discussed by Choksi and Ren, but both the result and the 
proof for this case are simpler than for the periodic cell: 

Theorem 3.3. Let (u, v) be a stationary point of F\ on R under constrained mass, with supp(u+w) 
connected and with a finite number of interfaces. Then 

1. Each pair of adjacent U-V type transitions is separated by the same amount; i.e. each U- or 
V -block is of the same width, with the exception of the two end blocks. 

2. The end blocks are half as wide as the internal blocks. 

The main tool in the proof of these theorems is the following lemma. 

Lemma 3.4 f [CR051 Lemma 4.1]). For any stationary point under constrained mass, the Poisson 
potential (f> has equal value at any two interfaces of the same type. 

The statements about the block sizes are deduced from this lemma, and from the fact that the 
potential cf> has prescribed second derivative on each block. 



3.3 Reduction to connected support 

We first need a technical result to rule out the possibility of an infinity of transitions. 

Lemma 3.5. Let (u, v) be a stationary point under constrained mass, let f2 be either R or Ti and 
let ui C ft be an open set such that v(uS) = {0}. Then lj contains at most two U-0 type transitions. 
A similar statement holds with u and v exchanged. 

Proof. On uj, </>" < 0; each U-0 or 0-U transition occurs at the same value of <j> (Lemma 13. 4p . say 
at <p = c G R. If the set {i£w: <f>{%) = c} has more than two elements, then by convexity, 

4>{x) = c for x G [xi,X2], 
</)(x) < c for x £ oj\ [xi, X2], 

for some x\ < X2 € to. On (xi, X2), therefore, 4>" = and thus u = 0. Therefore there are at most 
two transitions connecting U and 0, at x = x\ and at x = X2- □ 
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Corollary 3.6. If d uv > 0, then a stationary point under constrained mass has a finite number 
of transitions. 

Proof. By ©, at least two out of the three dij are strictly positive. If all three are positive, then 
the finiteness of Fi implies a bound on the number of interfaces. If one is zero, say d u $, then the 
lemma above states that the number of U-0 or 0-U transitions is no larger than the number of 
V-interfaces. Since the latter is bounded, the former is also. □ 

Theorem 3.7. Let N = 1. Let d uv > 0, and fix a mass M > 0. 

1. There exists a global minimiser under constrained mass M for which supp(u+v) is connected. 

2. This global minimiser is non-unique (apart from translation and mirroring) if and only if 

(a) the energy of this configuration is equal to the energy of another configuration (u,v) for 
which supp(u + v) is also connected, or 

(b) one of the following two conditions is satisfied: 

i. d u o = and there exists a global minimiser with an internal U-block or; 

ii. d v o = and there exists a global minimiser with an internal V-block. 

The non- uniqueness mentioned in condition I2al can manifest itself in multiple ways. Figure [2] 
shows how the optimal structure varies with mass: as the mass increases, the global minimiser 
progresses through structures with more and more layers. At the intersection points of the curves 
in the figure, indicated by a circle, structures belonging to different curves have the same value of 
the energy. Another possibility occurs when d u $ = d v o, since then u and v can be interchanged 
without changing the energy. The situation where two minimiscrs are both connected, have the 
same sequence of blocks (up to mirroring), but differ in the block widths, however, is ruled out by 
Theorem 13.31 

The fact that the global minimiser can be non-unique when, for example, d u a — can easily 
be recognised by an example. Suppose that there exists a global minimiser of the form UVUVU. 
Since the outer blocks of this structure are both U-blocks, Lemma T3.4I states that the value of <j> 
is the same at the two interfaces of U-0 type, and 4> is therefore symmetric around the middle of 
the structure. 

We now split the structure at the middle into two parts, and move the two parts apart. In 
doing so we create two new U-0 type transitions, which carry no energy penalty since we assumed 
d u o = 0. Since we split at the middle, where <j)' = 0, the new potential <fi can be constructed from 
the old one by translation of the parts, and the value of ||u — is also unchanged. 

Proof of Theorem \ S.7\ We defer the proof of existence of a global constrained minimiser to the end, 
and start by showing that existence of a global minimiser implies existence of a global connected 
minimiser. 

Suppose (u, v) £ K\ is a global minimiser such that M\supp(u + u) has at least three connected 
components. By Corollary 13.61 the support of u + v is bounded, and therefore we can take those 
three components to be (— oo,0), (2:1,2:2), and (2:3,00). The points 0, x\, x<i, and X3 therefore all 
are interfaces. 

Let <p be the associated potential; since u and v vanish on (x\, X2) and (0:3, 00), <fi is linear on 
(£1,2:2) an( A constant on (£3,00). Denote by <j>' 12 the value of <fi'(x) for x G [2:1,2:2]. 

For any < a < 2:2 — Xi, which we fix for the moment, we construct a new pair of functions u 
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and v with associated potential <fi as follows. Set 



(10) 



(11) 

x < Xl 

x\ < x < X3 — a (12) 
x > X3 — a 

Because <\> '(xi) = 4> (x\ + a) = 4>' 12 , the function <f> is continuously differentiable on R; and since 4> 
satisfies 0" = u — v on K, it is the Poisson potential associated with u and v. 

We now show that Fi(u,v) < Fi(u,v). As for the interfacial term in F\, if < a < X2 — x\, 
then the various transitions remain the same, only translated to different positions; therefore the 
interfacial term is unchanged. In the case a = X2 — X\, in comparison with (it, v) the two interfaces 
at x = X\ and x = x 2 have been joined to one interface, or have even annihilated each other; by 
the assumption §§§ this does not increase the interfacial term. 

For the second term of Fi we calculate 




We conclude that Fi(u,v) < Fi(u,v). Since (u,v) is a global minimiser, we conclude that 
Fi(u,v) = Fi(u,v) and thus that (u,v) is another global minimiser. Furthermore by Corol- 
lary supp(u + v) has a finite number of connected components and thus we can repeat this 
procedure until only one component remains. Therefore we have proved that if a global minimiser 
exists, then there (also) exists a global minimiser with connected support. 

Assume now that two global minimisers exist, one of which has connected support. The other 
global minimiser, let us call it (u, w), either has connected supp(u-fw) or disconnected supp(u + i'). 
In the former case we have proved part [2a] of the theorem; therefore we now assume the latter 
case, and show that this implies part I2bl 

Since (it, v) has disconnected supp(u + v), we can apply the construction above. For a given 
choice of a, we find another configuration (m, v) with energy equal or less than that of (it, v). 
Since (u,v) is a global minimiser, the energy of (it, v) is equal to that of (u,v) and thus the two 
inequalities encountered above are saturated. In particular, 

• The joining of the two interfaces surrounding a 0-block docs not reduce the energy; 

• The inequality (|13[) is saturated. 

The saturation of p^|) implies that <j)' 12 = 0, and therefore that (j>( x i) = 4'( x 2)- We now prove 
that these interfaces are of the same type, i.e. either both U-0 type or both V-0 type transitions. 

Suppose not, and to be concrete, suppose that the interface at x = x\ is a V-0 transition, 
and at 1 = 12 a 0-U transition. In this paragraph we will explicitly distinguish between mirrored 
interfaces of the same type, e.g. U-0 and 0-U. Since —0" = u — v and (j)'(xi) = 4>'{x2) = 4>i 2 = 0, 
there exists ay 2 > such that the next transition is at x 2 + y 2 and <j> decreases for x £ (x 2 ,x 2 +y2), 





u(x) 




X < Xl 


u{x) := < 


u(x - 


-a) 


xi < x < X3 








x > X3 — a 




f v{x) 




X < Xi 


v(x) := < 


v(x 4 


-a) 


Xi < X < X3 








x > X3 — a 










4>{x) := < 


<f>(x -1 


-a) 


— <f){x\ + a) + 



(f)(x 3 ) - (j)(xi +a) + 4>{xi) 
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implying that the next transition can not be a U-0 transition (which would require the same value 
for <p as at x — X2) but is a U-V transition, with a value of <p less than (j)(x2)- The same argument 
holds for the interface at x\\ the previous transition is at x\ — y\ for a y± > and is again a U-V 
transition, this time with a value of 4> larger than (f>(xi) = 4>(x2)- Since two U-V transitions have 
a different value of </>, the structure is not stationary, a contradiction. 

Since the interfaces at x\ and xi are of the same type, a non-changing interface energy implies 



that either d u o — or d V Q = 0, which is the first part of conditions 2(b)i and 2(b)ii Since the 
construction provides a global minimiser with an internal U-block (if d u o = 0) or an internal 
V-block (if d v o = 0), the second part of these conditions is also satisfied. 

We have now proved that existence of a disconnected global minimiser implies condition I2bl 
The opposite statement, that condition [2b] suffices for the existence of a disconnected global 
minimiser, follows from splitting any minimiser at a point x inside a U-block (supposing d u o = 0) 
such that 4>'(x) = 0. 

It remains to prove the existence of a global minimiser, and we now turn to this issue. Let 
(u n , v n ) be a minimising sequence. We first note that the translation arguments that we used above 
allow us to reduce an arbitrary minimising sequence to a minimising sequence whose elements each 
are connected. Therefore we may assume that the support of the sequence remains inside some 
large bounded set f2 C R, and does not approach the boundary of this set. 

Since both u n and v n are bounded in L°°(f2), there exist subsequences (that we again denote 
by u n and v n ) such that 

u n — ->> Woo and v n — u m in L°°(Q). 

Note that this convergence implies that J = J Vqo = M, since the constant 1 is an element 
of L 1 ^). Since L 2 {VL) C L 1 ^) we also have 

u n — Uoo and v„ — Voc in L 2 {VL). 

The functions u^, Voo, as the weak-* limits of u n , v n , take values in the interval [0, 1]. Thus if we 
replace K\ in (|2.2p . the definition of iq, by (note the change in value set) 

K\ := < (u, v) G (BV(fi)) 2 : u(x), v(x) G [0, 1] a.e., and uv = a.e., and u = v 



I Jn Jn 

then (uocWoo) G K\ and F\ is convex on L 2 (il). This implies that the subdifferential of F\ at 
(mocWoo) is non-empty, i.e. there exist pi,P2 G L 2 (Q) such that 

Fx(u n ,v n )> F(u 00 ,v 00 )+ / Pi(u n -u oa )+ / p 2 (vn-v 00 ). 

Jn Jn 

Weak convergence in L 2 (Q) now gives us lower semi-continuity with respect to this convergence: 

-Pi (uoo , Voo) < lim inf F 1 ( Un , v n ) . 

n — >oo 

It remains to prove that Moo and Woo are admissible, i.e. that they take values and 1 and that 
UooVoo = almost everywhere. In other words, we want to show that not only (uoo>^oo) G K\, 
but even (moo , Woo ) G K\ . 

By the assumption d uv > at least one of the coefficients c u and c v is strictly positive. Suppose 
that c u > 0; then the boundedness of J \u' n \ implies that the convergence of u n is strong in L 1 and 
pointwise almost everywhere |EG92( Theorem 5.2.4]. Therefore, for any ip G L°°(Q,), 

/ ipUooVoo = lim / ipu n v n = 0, 
Jn n ^°° Jn 

implying that Uoo^'oo = 0. Also the pointwise convergence gives 

Uoo G {0, 1} a.e. 
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If also c v > 0, then the same convergence holds for v^, and the proof is done. If instead cq > 0, 
then the same holds for u,^ + Voo, and again the proof is done. We continue under the assumption 
that Co = c v — 0. 

For the pair (uoo^oo) to be admissible, it is necessary that takes values in the boundary 
set {0, 1} only. This is a consequence of the lemma that we state below. □ 

Lemma 3.8. Let cq = c v = 0. If (u, v) minimises F± among all pairs (u, v) such that 

ii6 BV(R; {0, 1}) and v G BV(R; [0, 1]); 

• uv = a.e. in R; 

• Jr « = /r « = k u > 

then v(x) G {0, 1} for almost every 

Proof. Choose < rj < 1/2 and let w C R be the set of intermediate values 

uj = {x G K : G (77, 1 - 77)}. 
We need to prove that |w| = 0. Assume that |w| > and define a perturbation 

CO) = (00) - c)x w 0), 

where = U _^ is the Poisson potential associated with u — v, x^> is the characteristic function 
of the set lu, and c is a constant chosen to ensure that J £ = 0. Note that almost everywhere on oj 
the function </> is twice differentiable with <fi" > 77 > 0. 

Since the pair (u, v + e£) is admissible for e in a neighbourhood of zero, 

0= ^F 1 {u,v + eC) 



= 2 I Q4> = 2 I (0 - cY 

E=0 



so that <p is constant a.e. on u). As (f> is defined up to addition of constants we may choose = 
on uj. 

Since \uj\ > 0, we can choose xq G uj such that w has density 1 at xq and that <f> is twice 
differentiable at xq, with 4>"{xq) G (77, 1 — r?). Because of the density condition it is possible to find 
sequences a n G R, n G N, with the properties 

• a n — > as n — * 00; 

• For each n G N, xo ± a„ G w. 
Then 

0"Oo) = lim |a„|" 2 U(a:o - a n ) - 2(j>(x ) + (f>{x + a n )] = 0, 

n — ^cso 

a contradiction with <f>"(xa) > 77, and therefore with the assumption that uj has positive measure. 

□ 

3.4 Excursion: global minimisers on 

By very similar arguments one may prove the corresponding statement for functions on the 
torus Tl, thus extending the characterisation of CR05] to all global minimisers. 

Theorem 3.9. Let L > 0, d uv > 0, and fix a mass M > 0, with M < L/2. 

1. There exists a global minimiser (u,v) of F\ on the torus Tj, under constrained mass M for 
which supp(w + v) is connected. 

2. This global minimiser is non-unique (apart from translation and mirroring) if and only if 
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(a) the energy of this configuration is equal to the energy of another configuration (u,v) for 
which supp(w + v) is connected, or 

(b ) one of the following two conditions is satisfied: 

i. d u o = and there exists a global minimiser with an internal U-block or; 

ii. d V Q = and there exists a global minimiser with an internal V-block. 

Proof. The proof follows the same lines as in the case of R, Theorem 13.71 We will point out the 
differences between the two cases. 

Suppose (u, v) € Ki is a global minimiser such that Yl\ supp(w + v) has at least two connected 
components, which, by translating u and v, we can assume to be (x\, X2) C [0, L) and (#3, L) with 
X3 > X2- Let be the associated potential; since u and v vanish on (x^a^) and {x^,L), <fi is 
linear on these two intervals. Let 0'(x) = 0' 12 for x £ [xijXa] and cj>'(x) = 4>' 3L for x £ [x 3 ,L\. By 
possibly exchanging roles we can assume that |0' 12 | > \4>'^i\- 

Constructing for some < a < X2 — x\ the same translated functions u, v, and <p as given in 
(|T0l[T2|) . we have the analogous inequality 

+ <P' 3 l( L -X3 + a ) 
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■r 1 



(I 









/ 




J x\ -\-a 




/ 






J x\-\-a 


2 





<l <P' 2 +act>> 2 12 + ( f>>l L (L-x 3 ) (14) 

L 



Although (j) satisfies 4>" = u—v on (0, L), the function (f> can in general not be extended periodically, 
i.e. 0(0) 7^ </>(L). To correct this we define 

4>(x) :=0(x)-|(0(i)-0(O)), 

so that the function (f> solves cj>" = u ~ v on (0, L), is continuously diffcrcntiable on (0,L), and 
satisfies 0(0) = 4>(L). From 

4>'{L)-4>'{Q)= /V= f\v-u)=0, 
Jo Jo 

we conclude </>'(0) = 4>'(L), so that is the Poisson potential on associated with u and v. In 
addition, 

12 \{4>{L) 0(0)) f 0' + ^(0(L) - 0(O)) 2 







L 

,/2 



(I 



i(0(L)-0(O)) 2 



< /V- (15) 
Jo 

From these two inequalities it follows as in the proof of Theorem [XT] that F\ (u, v) < F\ (u, v), so 
that existence of any global minimiser again implies the existence of a connected global minimiser. 

We now turn to the discussion of the necessary and sufficient conditions for non-uniqueness. 
Again we use the fact that inequalities are saturated to deduce necessary conditions; in this case, 
however, there is an additional inequality in (|15p . The reasoning proceeds in two steps. 

Step 1: Take a < X2 — x\. When a < x-j — x\ no interfaces are created, annihilated or 
changed, and we only need to consider the inequalities in (fT4]) and (fT5|) . Since these are saturated, 
the following conditions hold: 
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1. |0' 12 | = \<j/ aL \, and 



2. 0(Z) = 0(0). 
We first calculate 

4>{L) - 0(0) = 0(x 3 ) - 0(xi +a) + 0(.xi) + cb' 3L (L -x 3 + a)- 0(0) 
= -<f>{xi + a) + (f>(xi) + acf>' 3L 
= a(03L - 0i 2 )- 

By condition [2] above we have 0' 12 = 3i > which is also compatible with condition [T] 

We now claim that 0' 12 = 4>' 3L = 0. Suppose not, say (for concreteness) 0' 12 = 4>' 3L > 0, then 
0(xi) < 0(x 2 ) and 0(2:3) < 4>(L). Since for a stationary point the potential has the same value 
at all U-0 type transitions and the same value for all V-0 type transitions, the two transitions 
at X\ and at x 2 are of different type, thus one is a U-0 type transition and the other a V-0 type 
transition. The same is true for x 3 and L (or 0); and the transitions at X\ and X3 are the same. 
Therefore 0(xi) = 0(2:3). 

For any fixed a in the interval (0, x 2 — xi), however, we have now constructed a second global 
minimiser (u, v) — and therefore a second stationary point — for which <fi(xi) ^= <p{x^, — a), since 
4>{xi) = 4>{x\) = 0(xi) and 

0(x 3 - o) = 0(x 3 - a) = 0(x 3 ) - 0(xi + a) + 0(x x ) 
= 0(xi) - a0' 12 
< 0(xi). 

Since the interfaces of (u, v) at x\ and X3— a are of the same type, this contradicts the stationarity of 
this second minimiser, and we conclude that 0' 12 = 3i = 0. Note that since the intervals (xi,x 2 ) 
and (X3, L) were chosen as arbitrary connected components of Tl \ supp(it + v), this implies that 
0' vanishes on the whole of \ supp(u + v). 

Step 2: Take a = X2 — x±. Non-uniqueness in this case implies that also the interfacial 
energy remains the same in the construction of (u, v). As in the case of R, the interfaces at x\ 
and X2 that are joined together in the construction of (u, v) are of the same type, i.e. either both 
U-0 type or both V-0 type transitions. The fact that is constant on 0-blocks is used in this 
argument. We conclude that either d u o — or d v o — 0, and that a connected global minimiser 
exists with at least one internal U-block (if d u o = 0) or at least one internal V-block (if d v o = 0). 
This proves the necessity of condition I2bl 

The sufficiency of condition [2b] follows by splitting one of the internal blocks, as in the case 
of R. Apart from the simplifying fact that the torus is bounded, the proof of existence of a global 
minimiser is identical to the case of R. □ 

Remark 3.10. Note that the proof of existence of a global minimiser generalises straightfor- 
wardly to the higher dimensional case of the torus , because the torus is bounded. On the 
unbounded domain M. N , N > 2, the above proof does not suffice. 

3.5 Explicit values and a lower bound 

We now focus again on functions on R. The results of the previous sections allow us to calculate 
global minima of the energy F\(u,v) as a function of the mass M = J u. Two important special 
cases are the monolayer and the bilayer. 

A monolayer consists of a single U- and a single V-block, of equal width m, where m is the 
mass of u or v, i.e. positioning the block around the origin for convenience, 

U(x) = X(- m ,0) and v ( x ) = X(0,m)! 
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where Xa 1S ^ nc characteristic function of the set A. We then find for the derivative of the Poisson 
potential 

for x < — m 
4>'{x) = < \x\ — m for — to < x < m 



for x > to, 



The total energy then becomes 



2 

monolayer of mass M = to: F\ = 2(cq + c u + c v ) H — m . 

3 

Note the definition of mass: a monolayer of mass M means that J u = J v = M, and therefore that 
the 'total' mass of the monolayer J(u + v) equals 2M. In this case the mass M of the monolayer 
equals the width to of each of the blocks. 

A bilayer consists of two monolayers joined back-to-back. It comes in two varieties, as UVU 
and as VUV. For a UVU bilayer of mass M = 2m, given by 

U{x) = X(-2m-m)u(m,2m) anc ^ v ( x ) = X(- m , m )j 
the derivative of the Poisson potential is 

for x < —2m 

—2m — x for —2m < x < —m 

4>'(x) = ^ x for —to < x < TO 

2m — x for to < x < 2m 

for x > 2m. 



The energy has the value 

4 
3 ? 

Similarly, n-layered structures consisting of n monolayers back-to-back, have energy 



4 , 

UVU bilayer of mass M = 2m : F\ = 2c Q + 4c u + 2c v + -to 
For a VUV bilayer the situation is of course analogous: 

4 , 

VUV bilayer of mass M = 2m : F i = 2c + 2c„ + Ac v + -to . 



vuvu. 


. V n- 


-monolayer 


with 


mass 


M = 


nm : 


Fi 


= 2c H 
= 2d v0 


- nc u + 
~t~ rid uv 


, „. 2n o 
(n + 2)c v + —m 3 

2n „ 

+ Y TO ' 


(16) 
(17) 


vuvu. 


.U n- 


monolayer 


with 


mass 


M = 


nm : 


Ft 


= d u o - 


- d v0 + 


2n , 
nd uv + —to , 


(18) 


uvuv. 


.U TV- 


monolayer 


with 


mass 


M = 


nm : 


Fi 


= 2d uQ 


~t~ nd uv 


2n q 


(19) 



Note that for a VUVU. . . V n-monolayer or UVUV. . . U n-monolayer the value of n is even, while 
for a VUVU. . . U rt-moiiolayer it is odd. Furthermore m is the U-mass in one monolayer, thus to 
is the width of the outer blocks, from which we see that the width of the inner blocks is 2to. By 
collecting these results wc find: 

Theorem 3.11. Let N = 1. For any structure of mass M , 

/9\ 1/3 / 

F\ > 2(c + min(c u , c)) + ( - J d^M. (20) 

In the limit of large mass, 

Jfmw{%ii :<-,.) 6*., /» = M} = (?) I/a e. (21) 
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Proof. If d uv = 0, then the first statement is easily checked and the second follows from the 
example of Section [2.21 We continue under the assumption that d uv > 0. 

Let (u,v) be a global minimiscr with connected supp(u + v), which exists according to The- 
orem 121 Note that for all three cases of structures (VUVU. . . V, UVUV. . . V, and UVUV. . . U) 
the intcrfacial terms are bounded from below by 2(co + min(c u , c„)), so that 

2 

F x > 2(c + min(c„, c v )) + nd uv + —rM 3 . 

3n z 

Minimising this with respect to n gives the desired lower bound. The particular value of n for 
which the lower bound is achieved, 

n (M) 3 := -— , 

will be useful below. 

To prove the second part of the theorem, we note that (|17H19|) imply the upper bound 
inf : («, «) e K U J u = m| < max{d„ , 4 } + inf | ^d uv + ^ (^j j . (22) 

Choosing the largest integer smaller or equal to uq(M) as particular value of n, 



n(M) 



'4 M 3 
3 duv 



= L^o(M)J, 



we have n (M) - 1 < re(Af) < n (M). In the limit A/ — > oo the quotient n(M)/M therefore 
converges to (4/3 d uv y^ 3 ; with this convergence the inequality (|2"2"]) implies (|2T|) . □ 

In Figure [3] the graphs depicting the energy per mass for VUVU. . . V configurations consisting 
of different numbers of monolayers are shown, for some specific parameter values. The lower bound 
from Theorem 13. Ill is indicated as well. 

Remark 3.12. Minimising F\JM from (fl6lIT9)) with respect to to, we find as minimising value 
for m, 

3 3(kid u0 + k 2 d v0 + nd uv ) 

m ° (n) := Tn ' 

where, depending on the configuration k\ = Q,k-2 = 2 (jTTJ) , ki =^2 = 1 p8[) or k\ = 2, k 2 = 
()19D . In all three cases we find that in the limit n — > oo, or equivalently (for to fixed) M — ► oo, 
the width of the inner blocks converges to 

2 lim TO (n) =6 1/3 di{ 3 . 

n — >oo 

Note that in |Mu l93] and |RW03aj it is found that one-dimensional minimisers for the functionals 
under consideration in those papers are periodic with period ~ (surface tension) 1 / 3 . (In these 
diffuse interface functionals the surface tension coefficients are given by integrating the square 
root of the potential.) 



4 Higher dimensions 

In this section we derive bounds on energy of minimisers in terms of the mass M. The first result, 
Theorem 14. 11 shows that the minimal energy has a lower bound that scales linearly in mass in the 
limit M — * oo. This is an extension of the lower bound (f2"0)) in one dimension, but with a smaller 
constant. 

A simple argument immediately gives an upper bound on the minimal energy at given mass: 
fixing any structure of unit mass, a candidate structure at mass M £ N can be obtained by 
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1.17 
1 



2 4 6 8 M 

Figure 3: Energy per unit mass for the one-dimensional case, according to the calculations in 
Section l3~ijl M is the total U-mass; for the parameters the values d u o = 1, d uv = 0.6 and d v0 = 0.4 
are chosen. All the graphs belong to a VUVU. . . V n-monolayer structure, where n/2 increases 
from 1 (left) to 20 (right) with step size 1. Also drawn are the (dashed) lower bound LB (|20p . and 

1 2 2. 

the asymptote 2 333^ ~ 1.17446. 

distributing M copies of the unit-mass structure over M. N . The energy of the resulting structure 
equals M times the energy of the unit-mass structure. This construction can be extended to 
non-integer mass M by spatially stretching a structure of integer mass close to M. In the limit 
M — ► 00 the resulting perturbation of the energy is small. 

In Sections l4.2l and l4.3l wc therefore provide tighter upper bounds, by constructing iV-dimensional 
structures out of near-optimal fc-dimcnsional ones, with k < N. 

4.1 Lower bound 

For this section we pick a function k G C™(R n ), non-negative and radially symmetric, such that 

/ k=1. 

In the following we will use the constant A jy, defined as 

A N :=i \e-w\dH N - 1 (w), 
J3N-1 

where S^" 1 is the (N — l)-dimcnsional unit sphere and e is any element of S' JV_1 . This definition 
is independent of the choice of e, because the integration is over all of S^" 1 . 

The central result is an interpolation inequality between the i31/-seminorm and H . In spirit, 
and in its application, it is similar to the Lemma 2.1 of [ChoOlj . The proof is different, however, 
and uses an argument of |KO02| , in combination with the characterisation of BV by |Dav02j . 



For e > we now define 



Note that L N k e = 1 for all e. 
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Theorem 4.1. Let d uv ^ 0. For all (u, v) G K\, 

f uKC^WWu-vW*^ Jf |Vu|Y , (23) 
where Ci(k,N) > is given by 

d(K,N) :=2iAj r [ 
The inequality \23\) also holds with u and v interchanged. Furthermore, 

F 1 (u,v)>C 2 (k,N) f u, (24) 

JR N 

where 

3 



Proof. If J RN u = the statements are trivially true. In what follows we assume J RN u > 0. 
First note that k £ * u G Hq (R n ) . From uv = it follows that w < 1 — u, so that 



(w — v)k 6 * u > / (2u — 1) ft e * u = 2 / mk £ * u 



Writing 



2 



/ uk £ * u = 2 / u(x)u(y)n e (x — y) dxdy 
Jr n Jr n Jr n 

{u(x) — u(y)) 2 K E (x — y) dxdy + 2 / u 2 
\u{x) — u(y)\n e {x — y) dxdy + 2 / u, 



R N JR N 



we have 



u< (u — v)n e *u+ / / — u(y)\n £ (x — y) dxdy. 

J R N JR N JR n 

The first term on the right-hand side is estimated by combining the definition of the if -1 -norm, 

(u — v)k £ * u < \\u — u||jj-i||Vk £ * u\\i,2, 
with the estimate (Young's inequality |Ada75( Theorem 4.30]) 

||V/s e *u|U» < \\u\\ L > / |VK £ | = ||«||ti / |Vk £ | =e' 1 \\u\\l 1 [ |Vk|. 
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For the second term we use a density argument as in |Dav02[ proof of Lemma 3] to find 
\u(x) - u(y)\n e (x - y) dxdy < 



< e 



Vu(ty+ (1 -t)x) 



(y - x) 



\y~x\ 



R N JO JR N 



Vu(x + th) ■ — 
| ft. | 



\y - x\ 



J — ^-K £ (h) dx dt dh 

£ 



- {x — y) dt dy dx 



VU(Z) ' \h\ 



- — ^-K £ (h) dz dh 

e 



\\Ju{z)-w\dH N - 1 {w)r N - 1 -K e {r) dzdr 
o Jr n Js n ~ 1 e 
r roo N 

eAnH*- 1 {S n ~ 1 ) / \Vu{z)\dz / —K £ {r)dr 
Jr n Jo £ 

£A N / \y\n{y)dy I \Vu(z)\dz. 



The first equality follows after substituting y = x + h, while the substitution x — z — th leads to 
the second equality. 

Collecting the parts we find the estimate 



U < £ 1 \\U~V\\ H -1 



|V«|+eC (K,iV) / |Vu|, 



where 



C (k,N) :=A n / \y\K(y)dy. 

IR N 



Minimising the right hand side with respect to e we find 

/ 

Jr n 



u<2 







C (k,N) f \Vk\\ 


/ u ) \\ u -v\\h-i 


Jr n ' 


Ur n ) 



IVul 



Dividing both sides by (J mN u) 4 and then raising them both to the power 4/3 gives the first 
statement of the theorem. Since we have used no property that distinguishes u from v, we can 
apply the same argument with u and v interchanged. 

To prove the inequality |24|) , we remark that from (|23|) and Young's Inequality we obtain, for 
any a, (3 > 0, 



Cf7»<f/|V«l + 3>- 



By choosing 



-}J Z and (3 



a/3 



and then adding the two inequalities with weights a 2 and j3 2 respectively, estimate ([24)) follows. □ 
Note from the proof above that estimate d24|) is not sharp if L. N u > 0. 



Remark 4.2. Inequality ([24]) does not hold in the case where d uv = 0. The same sequence 
(u n ,v n ) that was introduced in Section 12.21 demonstrates this fact, since Fi(u n ,v n ) — > 2 while 
/ u n — > 00. 
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4.2 Upper bound 

We next show that the one-dimensional upper bound (f2"Tj) (or ([3])) also holds in higher dimensions, 
as a consequence of the more general statement below. Theorem 14.31 formalises the intuitive 
idea that extending a one-dimensional minimiser in the other directions, and then cutting off the 
resulting planar structure at some large distance, should result in an N- dimensional structure 
whose energy-to-mass ratio is close to that of the original one-dimensional structure. We formulate 
the result for fc-dimensional structures that are embedded in N dimensions. 

Let 1 < k < N — 1, and let us write Ki >k for the admissible set K\ on R fe . Let (u,v) be 

• any element of K\^ k , when k > 3; or 

• any element of with J Rk x(u(x) — v(x)) dx = 0, when k £ 1, 2. 

(We explain this restriction in Remark 14.5ft . Split vectors x € WL N into two parts, x — (£,??) £ 
R fe x K JV-fc ; and dcfine a cutoff f unction Xa . ^N-k [ Qj ]_] by 

Xa(v) ■= x(\v\ - a ). 

where \ : R — ► [0, 1] is fixed, smooth, and satisfies x( x ) = 1 for a; < 0, x( x ) = f° r x>l. We will 
compare the energy values of the k- dimensional structure (m, v) with those of the N- dimensional 
structure 

(u,v)(x) := (u,v)(0Xa(v)- (25) 

Note that this (it, v) is an element of -K^jv, the admissible set K\ on M. N . 

A note on notation: LOd will denote the d-dimensional Lebesgue measure of the d-dimcnsional 
unit ball. 



(26) 



Theorem 4.3. Fix (u,v) as given above. Then, for (u,v) as defined in &25\) . 

F 1 (u,v) Fi(w,u) . nn/ , 
P = — t. — — \-0(l/a) as a — ► oo. 

JR N U JR k U 

Proof. We first estimate the interfacial terms as follows: 

|Vu|=/ / \Vu(0\x a (ri)dr)dt+ [ f u(0\^X a (v)\ <M£ 

JR k JR N - k JR k JR N ~ k 

< Lo N . k {a + l) N ~ k f \Vu\ + (N- k)u N - k (a + l^-^llx'Uc 

jR k 

> w N ^ k a N - k [ \W\, 

< jR k 

and therefore 

|Vu| = (1 + 0{l/a))oj N - k a N ~ k / |Vu| as a -> oo. 

R v JR k 

Similarly, 

f \Vv\ = {l + 0{l/a))u N ^ k a N - k [ \Vv\ and (27) 

JR N JR k 

|V(u + «)| = (1 + 0(l/a)V JV _ fc a iV - i: / |V(u + u)| (28) 

R N JR k 

The estimate of the _1 -norm is formulated in Theorem l4.4l The result now follows by combining 
the estimates (|2l)H2l?|) and remarking that the mass of (u, v) is given by 

/ u = f u(0Xa(v) = (1 + 0(\/a))u N . k a N - k f u(£). 
Jr n Jr n JR k 

□ 
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Theorem 4.4. Under the conditions above there exists a constant C = C(k,N) such that for all 
a > 0, 



<Ca N - k - 1 I [|V^| 2 +/]. (29) 



\ u ~ w lii?-i(R«) ~ u N _ k a N k \\u - v\\%- HRh) 
Here <f> is a k- dimensional Poisson potential associated with (u,v). 

Remark 4.5. The restriction of vanishing first moments for k = 1,2 follows directly from 

2 

the requirement that J Rk <f> can be chosen finite in (f2"9"|) . Since the integral of u — v vanishes the 
potential <f> := G*(u — v) decays to zero at least as fast as |£| 1_fe , as can De seen horn the multipolc 
expansion of <f> (see [HKDS93] ). For dimensions k > 3 it follows that J 4> is finite; but for k = 1,2 
a higher decay rate is necessary, which we provide by requiring an additional vanishing moment. 
The case k = 1 is special: the vanishing of the zero and first moments implies that <fi :— G*(u—v) 
is zero in a neighbourhood of infinity. 

Proof of Theorem \4-4\ The Poisson potential <j> associated with (u,v) satisfies 

-Acf>(x) - u(x) - v(x) = («(£) - v(0)Xa(v) for x = (f , rf) G R N . 
Similarly, the fc-dimcnsional potential <fi associated with (u, v) satisfies 

We write for the part of the gradient that operates on £, that is (d xi , d X2 , . . . , d Xk , 0, . . . , 0), 
and we use a similar notation for the other part of the gradient V„ and the partial Laplacians 
and A,,. Remarking that 

Xa(»?)V0(x) • V^(£) dx = [ Xafa)V^(ar) • V^(£) dx = - [ X a(v)<t>(x)A e ^) dx 

Xa(ri)<t>(x)(u(0 - v(0) dx = - [ <j){x)A4>{x) dx 

JVL N 

|V0(x)| 2 dx, 



we calculate 



/ |V0- X aV0| 2 -/ |V0| 2 -2 / XaV0V0+ / Xa|V0| 2 

Jr n Jr n Jr n Jr n 

|V0| 2 +/ xlm 2 - (30) 



One inequality relating the two norms can be deduced directly: 

\\u-vf H - Hm) = ( |V0| 2 < /" x 2 a \VW <u N . k (a+l) N - k \\u-v\\ 2 H . im 

For the opposite inequality we set 

t/j(x) := <f>(x) - ~${0Xa(v), 

and rewrite 

|V0- X aV0| 2 - / |V^| 2 + 2/ 0V^V, ;X a+ / 0Vxa| 2 

|W| 2 -2 / V0A I)Xq + / 2 |VXa| 2 

JR™ Jl» 
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Since 



-Alp = -A(j> + Xa\4> + 4>Ar)Xa = Xa(u - v) - \a{u - v) + 0A r ,Xa = (/>\Xa, 

the function i\) is the global minimiscr of /, which is convex as functional on ip. Therefore, setting 
M*) -^(e)|Vx (r?)| 2 , 



|V0| 2 |V Xa | 4 + 40 2 |D a X « ■ V X a| 2 - 20 2 |VXa| 2 Axa + 0lVxa 12 



<C( X )(a + l) 



N-k-l 



|V0| 2 



where the constant in the last line depends on x but can be chosen independent of a. 
Combining this estimate with (|30[) provides us with the opposite inequality, 



uj^]^kO N k \W — V " 2 



H -, m < / X 2 |V0| 2 

— Il u — ' 1, IIh-1(R«) + 

+ C( X )(a + l) N - k - 1 



(31) 



Summarising (j3"Tj) and (|30|) as 



iv-fc ^w-fc 



N-k\ 



\\ u - W|| H -i(K«) - ^JV-feO ||U - UHjj-ifRk) 

we find the statement of the lemma. 



< uJN-k((a + 1) 

>-C( X )(a + l) N - k -iJ Rk |W| 2 +? 



a ) \\u — v 

.-r,o - 



4.3 Examples with prescribed morphology 



□ 



Theorem 14.31 has the following consequence: when comparing energy-per-unit-mass of structures 
in dimension N, we can include the energy-per-unit-mass of structures in dimension k < N, up to 
a correction term that decays to zero in the limit of large mass. We now use this tool to investigate 
the energy values of various fixed-geometry structures. 



• A one-dimensional, lamellar, structure. The optimal energy-per-unit-mass is , 
(Theorem 13. Ill achieved in the limit of large mass). 



9 A 1/3 ,2/3 



• A micelle in N dimensions, i.e. a spherical particle described by 

w m (a;) := 

v m (x) := 

The equal-mass criterion implies that R2 = 2 1 / N R\, and by optimising with respect to the 



1 if < \x\ < Ri, 

otherwise, 

1 if R t < \x\ < R 2 , 
otherwise. 



remaining parameter R\ we find that the optimal energy-per- unit-mass is (Theorem IA. 1 (1 

1 

V 



3 (d uv + d v0 V2j 3 ( log 2 



for N 



and 



for N > 3. 



(32) 



These optimal values are attained at finite mass. 

In both cases the micelle energies are larger than (f) 1 ^ 3 dtl? , even when d V Q = 0, implying 
that for large M lamellar structures have lower energy per unit mass than micelles. 
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• A /c-dimcnsional micelle embedded in iV-dimcnsional space, similarly to the case of lamellar 
structures in TV dimensions. A two-dimensional micelle thus becomes a cylinderical structure 
in three dimensions. 

The energy per unit of mass of such a structure will be lower than that of an iV-dimensional 
micelle, since (|52"j) is a strictly increasing function of N, but larger than that of a lamellar 
structure for large mass, by the conclusion of the previous point. 

• A monolayer in the shape of a spherical shell as in Theorem IA.ll Here the optimal energy 
per unit mass can be found (in the limit of large radius R) by minimising (|33[) with respect 
to M : 

1/3 

(d u0 + d m + d v of /3 + (N - 1) (dyQ - d u0 ) 

+ V3 (n - 1) (-^^Ko + d v0 ) + 4») ir 2 + o(ir 3 ). 

Note that this value approaches for R — > oo the optimal one-dimensional value when d u o = 
d V Q = 0. (Although such a choice is ruled out by ©, one may calculate the value of the 
energy in this case nonetheless.) In this case the limit value is approached from above. 

Alternatively, if either d u o or d v o is non-zero, then the limit value is larger than that of the 
optimal lamellar structures. 

Among this list, therefore, the structures with lowest energy per unit mass are the lamellar 
structures. It seems natural to conjecture that global minimisers also resemble cut-off lamellar 
structures, and have comparable energy per unit mass. On the other hand, the results of the 
companion paper [GP07j show that bilayer structures are unstable in a part of parameter space, 
and similarly Ren and Wei showed that for the pure diblock copolymer model (u + v = 1) 'wriggled 
lamellar' solutions may have lower energy than straight ones |RW05] . Determining the morphology 
of large-mass global minimisers is therefore very much an open question. 

5 Discussion and conclusions 

The results discussed in this paper provide an initial view on the properties of the energy |T]) , and 
consequently on mixtures of block copolymers with homopolymers. The sharp-limit version of the 
more classical smooth-interface energy provides a useful simplification and provides us with tools 
that would otherwise be unavailable. 

In one dimension we continued on the work of Choksi and Ren and gave a complete character- 
isation of the structure of one-dimensional minimisers, both on R and on a periodic cell. 

In the multi-dimensional case we have proved upper and lower bounds for the energy of min- 
imisers. These bounds both scale linearly with mass, but have different constants. The upper 
bound is derived from the one-dimensional minimisers, thanks to the cut-off estimate of Theo- 
rem |4j3j the results of the companion paper on the stability of mono- and bilayers |GP07| suggest 
that for some parameter values this upper bound can be exact, while for others it is not. Similarly, 
the lower bound proved in Section [4. II has the right scaling in terms of mass, but the constant is 
not sharp. 

The sharpness of the estimates is especially relevant in relation to the issue of optimal mor- 
phology. A precise estimate of the energy level of energy minimisers may exclude large classes of 
morphologies and thus limit the possible morphology of energy minimisers. 

Since we lack such a sharp estimate the question of the preferred morphology in multiple 
dimensions is still completely open. Part of this question is the behaviour of the morphology near 
the copolymer-homopolymer interface. For instance, if the preferred morphology is lamellar, does 
the lamellar orientation show a preference to be orthogonal, parallel, or otherwise aligned with 
respect to the interface? The experimental observations of for instance |KHH94( ION97( [ZJM05 
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show both orthogonal and parallel alignments. Other issues are those of the penalty incurred 
by certain macrodomain morphologies and defects (such as chevron or loop morphologies [AM04 . 
Figures 11 and 14]). 

The large-mass limit for the functional F\ is equivalent to a singular-limit process at fixed 
mass for the functional F e ([6]). As discussed in Section IT~5l the results of [PR06] suggest that for 
certain values of the Cj — to be precise, for those values for which bilayer structures are stable — the 
functional F £ may display similar, partially localised behaviour. 

The results from this paper have some interesting physical implications. Theorem 14.31 tells 
us that extended one-dimensional minimisers, i.e. layered structures, will have a relatively low 
energy (although the question whether or not these are minimisers is still open). Theorem 13.31 
and Remark 13.121 show that these layers all have the same width. One can think of lamellar 
configurations like this as having all polymer molecules aligned in straight rows next to each other. 
Structures like 0U0, which are not to be expected on physical grounds, are a priori not forbidden 
in our model, but such configurations are not stationary points, as is shown in Section [3.31 

Depending on the surface tension coefficients very different structures can appear. As remarked 
in Section 12.21 the role of d uv is a special one. If d uv = 0, there is no repulsion between the U- 
and V-phases, but there is attraction, due to the if -1 -norm. Complete mixing of both phases 
will occur. If one of the other surface tension coefficients is zero instead, say d u o = 0, then UVU 
bilayers can be joined together without extra cost, and vice versa UVUVU structures can be split 
through the middle without increasing the energy. Physically this happens if the U- and 0-phases 
do not repel each other. The simplest case one can think of is if both phases consist of the same 
material. If both d u o > and d v0 > 0, it will always be energetically favourable to join different 
layers together, because doing so decreases the length of the energetically costly interfaces. 

A Spherically symmetric configurations 

In this appendix we will compute the energy F± of spherically symmetric monolayers and bilayers. 
In |ON98| the energy for a spherically symmetric bilayer in two and three dimensions is computed. 
We will give the energy in any dimension N. 



A spherically symmetric monolayer with inner U-band in N dimensions consists of a spherical 
layer of U between radial distances Rq and R\ and a spherical layer of V between radial distances 
R\ and i?2- An example for N = 2 is drawn in Figure l4al Similarly, a spherically symmetric UVU 
bilayer is a spherical layer of V between radial distances R\ and R2, flanked by two spherical 
layers of U, between radial distances Rq and R\ , and R2 and R3 respectively. A two-dimensional 




(a) A spherically symmetric monolayer in 
two dimensions 



(b) A spherically symmetric bilayer in two 
dimensions 



Figure 4: 
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example is shown in Figure I4bl Monolayers with inner V-band or VUV bilayers are constructed 
by interchanging U and V. 

We will first compute the energy for monolayers. In Theorem IA.1I we give the expansion in 
terms of the curvature k of the energy per mass, for small k. The exact expressions for F\ can be 
found in the proof of the theorem, in ([35)) for N = 2 and in ([36)) for N > 3. 

The expansion in terms of small curvature is obtained by linearising these exact energy expres- 
sions around ll R = oo" . To this end we introduce for the monolayer the curvature k, total U-mass 
M, and mass per (hyper-)surface area m: 

K:=iJf\ M := u N (R? - Bg), 
M Mk n ~ x 



We then get 



Rq = k _1 \/l — Nmn, 



R 2 = k" 1 \/l + Nmn. 

Theorem A.l. Let (u,v) € K\ be a spherically symmetric monolayer with inner U-band. Fix the 
mass per surface area m > 0, then for all N > 2: 

F\ 2 

— (u, v) = m" 1 (d u0 + d uv + d v0 ) + -m 2 + (N - 1) (d v0 - d u0 ) k 

+ (N- l)m (~(d u0 + d v0 ) + -^(3N - 2)m 3 ^) k 2 + 0(k 3 ), (33) 

if k io. 

Note that there are two configurations for a monolayer, depending on whether the U-phase 
is on the inside or the outside. The theorem above states the case where the U-phase is on the 
inside. The other case is found by interchanging d u o and d v o- 

Proof of Theorem \A.1\ The proof consists of three steps. First we compute F\ in terms of the radii 
Ri, then we rewrite it as a function of k, M and m. Finally the expansion is found by computing 
the first terms of the Taylor series of these expressions for k -C 1. 

The interfacial terms are computed in a straightforward manner. For the iJ _1 -norm we need 
to compute the Poisson potential, which depends only on the radius r because of the spherical 
symmetry and which we will denote by <fi(r). The Poisson equation in spherical coordinates is 

_ r -JV+i (yv-i^r))' =u-v for r > 0, 

#o) = ^(o) = o. [6 ' 

The solutions to this equation for N = 2 and N > 3 are different, and we treat these two cases 
separately. First we solve for N = 2 for the four different regions and we match the solutions 
under the condition that <fi G C 1 (R Ar ). This gives 



(r) 



re(0,Ro), 
-y + lR 2 logr + \R 2 - ii?ilogi? r e (R ,Ri), 

i- 2 J- | (jRg - 2R 2 ) \ogr + R\ logi?i - \R 2 \ogR a - \R\ + \R 2 re (R u R 2 ), 



(2R 2 -R 2 )(l- log(2i? 2 - R 2 )) + R\ log R X — \R^ log Rq — + jR% r> R 



2- 



Note that (j) is constant on [R2, 00): the solution can not have a term proportional to logr on this 
interval, since 4> € W 1 ^ 2 (R N ) and 



/>oo />oo / ~t \ 2 

/ \d r \ogr\ 2 rdr = / (-) r dr = 00 

JR 2 JRi \ r J 
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This means that 0'(i?o) = 4>" '{R2) = 0. We now compute the norm via 

rRi />-R2 

<f)(r)r dr — I 

I R JR.! 

For N = 2, we then find 



w ~ w llff-i(fi) = 2?r I / <t>( r ) rdr - I 4>{r)rdr 



v) = Roduo + Rid uv + R%d v0 

Z7T 

- Jitf + ii?gi?2 _ li?* l 0g i? _ i?2 ^2 _ jjg) lQgRi 

+ ±{2R 2 -R 2 ) 2 \og(2R 2 -R 2 ), (35) 

where the radii are related by R 2 , — R\ = R\ — R 2 ,. 
Analogously solving for N > 3 we find 

ifre(0,i? o ), 

-1 nW '-N+2 _ 1 _3 1 1 p2 tc _ c c p p \ 

N(N-2) n Q ' 2JV' 2(W-2) rl 11 ' fc K^Ot^Dt 

} = i iv^b) « - 2 <) ^ + W^ 2 + spby K - 2iZ?) if r G (i? l5 R 2 ), 
5^ ((2i?f - <) * - 2i?f + i$) if r > i? 2 , 

and compute the norm. This leads to 



— (u, v) — Rq 1 d u o + R± 1 d uv + R2 1 d V Q 



«+ 2 - <+ 2 ) + ^r^Ri « - <) , (36) 



N 2 — 4 y 2 ' N(N -2)' 

where the radii are related by R2 — Ri = R\ — Rq ■ 

Rewriting our results in terms of k, M and m gives, for N = 2, 

p. 1 

— (u, v) = Ar 1 (yi - 2mnd u0 + d uv + \/T+ 2mnd v0 ) - -n~ 2 

+ i f im _1 «T 3 + kT 2 + log(l + 2mn) 

- i ^im _1 K -3 - k 2 + tokT 1 ^ log(l - 2mn). 
For the monolayer with N > 3 we get 



■^( u i v ) = m 1 ((f _ Nmn^N 1 d uQ + d uv + (1 + JVm«) ^d„o) 



1 1 /, , JV + 2 , , JV + 2 \ 

j^— -m- 1 (Jl - Nm K )— - (1 + A^m K ) — J k " + 



.-3 1 * ^-2 



We now expand in terms of k <C 1 to get the result. □ 



Next we turn to the bilayer. Here we follow the same route as before. Theorem IA.2I states 
the expansion in small curvature; the exact expressions for F\ can be found in the proof in (1381) 
for N = 2 and in (J35J) for N > 3. Define R > via R N = \ + R%), then for the bilayer we 
introduce curvature k, total U-mass M and mass per (hyper-)surface area m for N > 1 as follows: 

K:=ir\ M := u N {R% - i?f ), 
m := — . 
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Then 



Ro 



^{l-Nmn), Ri = K 



1 \ 
—Nmn , 

2 r 



R 2 




R 3 



k" 1 V(l + Nm K ). 



Note that here we have chosen the radii such that the inner and outer U-band have equal mass: 
2 



R\ 



Rf 



i?Q. This is not the optimal choice, in the sense that a small change in the relative 



thicknesses of the inner and outer monolayers might improve the energy slightly. We expect this 
to be a small effect, however. 

Theorem A. 2. Let (u,v) € K\ be a spherically symmetric UVU bilayer. Fix the mass per surface 
area m > 0, then for all N > 2: 



Fi{u,v) 
M 



2m- 1 (d u 



>)+a m2 
o 
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^(3A-2)m 3 j K 2 +0( K 4 ), (37) 



An analogous result and proof corresponding to the VUV bilayer is constructed by replacing 
d U Q by d v0 . 

Proof of Theorem ] A. 21 As in the proof of theorcm lA.ll wc follow three steps. First we compute F\ 
in terms of the radii Ri. The resulting expression we rewrite in terms of k, M and m and finally 
we find the expansion in terms of k <C 1 by calculating the first terms of a Taylor series. 

The main problem in the first step consists of deriving the Poisson potential that solves (f34]) . 
For N = 2 we find 
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if r > i? 3 . 



For JV > 3 we have 
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We then proceed in the same way as for Thcorcm lA.il to find, for N = 2. 



if re (0,.Ro), 
if re (Ro,Ri), 
if r e (i?i,i? 2 ), 
if r e (i? 2 ,i? 3 ), 
if r > i? 3 . 
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where R 2 ~R 2 = ±(i? 2 - R{) = R\ — R%. 
For a bilayer with N > 3 we have 

- -Fi (u, v) = (iC" 1 + R^duo + (R?- 1 + + « +2 - < +2 ) 



Nw N ^ ' ' v u d ' " u v 1 ' ' 2iV(iV + 2) 

'-4<i? 2 + 4<i? 2 - 8<i? 2 + <+ 2 + 4<+ 2 + 4<+ 2 - <+ 2 ) , (39) 



2N(N-2) 



where i?f - R? = \(R$ - iif ) = i?f - i?^. 

Rewriting these results in terms of k, M and m gives, for N = 2 
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For > 3 we have 
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The result (|3T|) now follows from expanding in k< 1. □ 
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